— USSCI  ASS.IF.1.T& 

'"t  A'-''*  >'  A »sfis  op  this  t**r i /up,.  » *»*»?*-.<) 

I REPORT  DOCUMENTATION  PAGE 


A!)  INN  I klif  i ION*? 

ni  rnwj  { 'ivpt  j *jNr,  i*oru 


;*  GOVI  ACCESSION  NO  j 3 J<CClR»*  M?*S  CATALOG  NUMBER 


cr /r**t /nr  „ 

'Ll  i*)^S»6'  l»> 


A CAICULAT lOiiAL  MODEL  FOR  HIGH  ALTITUDE  FHP 


ZLOg 


I % type  op  RL^oflT  r,  prnioo  covered 


*6  PERFORMING  ORG  REPORT  NUMBER 


Louis  W,  Seiler,  Jr. 

Capt  USAF 

h f7rVoRM,NG  ORGANIZATION  name  ano  address 

Air  Force  Institute  of  Technolony  (AFIT-EN'P) 
Wrirht-Patterson  Air  Force  Base,  Ohio  ^5^33 

7i  CONTROLLING  OPPlCE  NAME  ANO  ADDRESS 


| 8 CONTRACT  OR  GRANT  NUMBERS) 


| 12  REPORT  DATE 

I March  1^75 


Air  Force  Weapons  Laboratory  (DYT)  March  1Q75 

Kirtland  AFB,  New  Mexico  87115  ^“number  op  pages 

i C*o 

* MCNJTCRlNG  AGENCY  NAME  ft  AOORfSS.  i ditlerenl  (tom  Controlling  0//icc,  i 15  SECURITY  CLASS  (ol  thit  repott) 


Air  Force  Institute  of  Technology  (AFIT-ENP) 
Wr ieht-Patterson  Air  Force  Base,  Ohio  A5**33 


UNCLASSIFIED 


] 1ft  OlSTRjBuTiCN  statement  rot  this  Report) 


Approved  for  public  release;  distribution  unlimited. 


| 17.  DISTRIBUTION  STATEMENT  (of  the  obst-act  wtwfiir  Clock  20,  il  dltterent  (tom  Rrpoit) 


I 16  SUPPLEMENTARY  NOTES 


Approved  for  nubile  release, IAW  AFR 
, ‘ 190-17 

Jeriy  C.  Hix  Captain,  USA? 

Director  of  Information 


KEY  *0R0$  rContmoo  or,  revet  so  it  Jo  ,(  nccer'<**>  wd  /deo*»/>  k\  Mock 

EMP 

High  Altitude  EMP 
Electromagnetic  Pulse 
EMP  Model 


(•produced  by 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 


25.  ABSTRACT  (Continue  on  rovers*  tide  Hn«tiN,j  and  Identity  by  btoc*  nuvber) 

An  electromagnetic  pulse  (EMP)  model  is  developed  which  allows  a quick 
computation  of  the  tine  development  of  the  electric  fields  oenerated  hy  a high 
altitude  nuclear  burst.  The  model  is  based  on  the  Karzas-Latter  hinh  frequency 
approximation  for  high  altitude  EMP,  which  describes  fields  generated  by  fomoton 
electrons  interactlno  with  the  earth's  rannetic  field.  The  proper  choice  of  a 
eam-ia  time  output  function,  which  can  be  integrated  in  closed  form,  and  a small 
ancle  approximation,  made  in  the  expressions  for  the  Compton  currents  and  air 
conductivity,  eliminate  the  time  consuming  numerical  inteerntlons  usually 

DD  , j2n”j  1473  eoit.™.  or  !NOv«s  is  obsolete  * UNCLASSIFIED 

' s E CU ° I T Y CLASSIFICATION  OF  This  PAGE  (Hhen  Duo  JTnfrrcd; 


'ecessary  in  IMP  models  to  connute  the  Conotcn  currents  and  air  conductivity. 
"His  results  in  a ccnsiviornblc  savinns  in  co- nut  it  ion  ti'*’**.  The  rode  I is 
""esenteJ  in  a nt*r  ihicb  is  sinr»le  to  us^  hut  still  nMo“S  the  variation  of 
t ’2  rt'VOi-  th  oret ica I nar.neters  in  tlie  nrohier.  & simplifies  model  of  clectror 
collision  frequency  as  a function  of  electric  fieid  strennth  is  pivcn  which 
''nobles  the  ncdel  to  predict  accurate  results  for  nuclear  weaoon  n«w?  yields 
uo  to  at  least  100  Kt.  The  results  predicted  by  this  EMP  model  compare  to 
-i thin  S.Sr'  with  results  fron  the  Air  Force  Weapons  laboratory  CHEMP  computer 
cc^e.  The  computation  time  us  inn  the  presented  rode  I on  a COC  Ef.00  comnuter  Is 
tvnically  5 see  o.-  less  for  a 5 shake  computation  ner tod  in  steos  of  .1  shake. 
Tne  model  presented  should  be  useful  for  both  classroom  instruction  and 
"uclear  vulnerability/survivability  studies  and  analysis  problems. 


if IWCLASSlF  Igf) _ 

SECUmtV  CLASStFlCATtON  OF  ThIS  PAGCnThcn  Er,ff4) 


A calculational  model  for  high  altitude  emp 


GEP/PH/75-13 


Louis  W.  Seller,  Jr. 
Captain  USAF 


D D C 

j\P'~'r  1nn  rjrj 

j -;  m « H15 

IteSiaiuTS 

B ^ 


Approved  for  public  release;  distribution  unlimited. 


1 


A CALCULATIONS.  MODEL  FOR  HIGH  ALTITUDE  EMP 


Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 

In  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Haster  of  Science 

by 


Louis  W,  Seiler,  Jr.,  B.S. 
Captain  USAF 


Graduate  Engineering  Physics 
March  1975 


Approved  for  public  release;  distribution  unlimited. 


ti» C 


GEP/PII/7S-I3 


Preface 


I would  like  to  express  my  appreciation  to  my  advisor,  Major  Carl 
T.  Case,  for  his  patient  guidance  and  understanding  during  the  work  on 
this  thesis.  His  advice  was  a major  contribution  towards  the  comple- 
tion of  this  work. 

I am  also  gratefully  indebted  to  Captain  Leon  A.  Wittwer  of  the 
Air  Force  Weapons  laboratory  for  his  enlightening  insights  Into  the 
finer  aspects  of  the  high  altttude  EHP  problem. 

I would  also  like  to  thank  Captain  Terry  C.  Chapman  for  his  time 
in  discussing  various  concepts  and  ideas  with  me. 

A large  amount  of  my  thanks  must  go  to  my  wife,  Joan.  Without 
her  patient  understanding,  moral  support,  and  encouragement,  this  thesis 
could  not  have  been  completed. 


Louis  W.  Seiler,  Jr. 


"O. 


GEP/PH/75-13 


Contents 


Page 


II. 


III.  Model  Considerations 


Preface II 

List  of  Figures  Iv 

Abstract  v 

I.  Introduction  1 

Theory  4 


Overview  .........  4 

Particle  Densities  4 

Currents  and  Conductivity  8 

Field  Equations 12 


IS 


General IS 

Tine  Dependence  of  Gamma  Rays  IS 

Energies 17 

Compton  Electron  Lifetime  22 

Small  Angle  Approximation  22 

Electron  Collision  Frequency  23 

Preionization  25 

Calculation  Method  26 


' VYWt- 


/ 

\ 


JL 


i 


4 


G£P/f H/75-13 


List  of  Ftguies 


Figure 

1 Model  Geometry  ................... 

2 Typical  Lone  Gamma  Output  Pulse  

3 Typical  Narrow  Gamma  Output  Pulse  

k Peak  Electric  Field  Values  for  100  Km  Burst  Heights 

5 Electric  Field  Development  for  .25  Kt  Ganna  Burst  . 

6 Peak  Electric  Field  as  a Function  of  Burst  Height  . 

7 Peak  Electric  Fields  for  Two  Gamma  Pulse  Shapes  . . 

8 Effect  of  a .03  Kt  Gamma  Yield  Precursor  Burst  . . . 

9 Effect  of  Geomagnetic  Field  Strength  on  Peak  Fields 


Pane 

. 5 
. 18 
. 19 
. 30 
. 32 
. 33 
. 3*i 
. 36 
. 38 


GEP/PH/75-13 


l 


| 

l 


I 


l 


Abstract 

An  electromagnetic  pulse  (EMP)  model  Is  developed  which  allows  a 
quick  computation  of  the  time  development  of  the  electric  fields  qene- 
rated  by  a high  altitude  nuclear  burst.  The  model  Is  based  on  the 
Karzas-latter  high  frequency  approximation  for  high  altitude  EMP,  which 
describes  fields  generated  by  Compton  electrons  interacting  with  the 
earth's  magnetic  field.  The  proper  choice  of  a gamma  time  output  func* 
tion,  which  can  be  integrated  in  closed  form,  and  a small  angle  approxi- 
mation, made  in  the  expressions  for  the  Compton  currents  and  air  conduc- 
tivity, eliminate  the  time  consuming  numerical  Integrations  usually 
necessary  in  EMP  models  to  compute  the  Compton  currents  and  air  con- 
ductivity. This  results  In  a considerable  savings  In  computation  time. 
The  mode'  Is  presented  in  a manner  which  is  simple  to  use  but  still 
allows  the  vari  ’■'on  of  the  major  theoretical  parameters  In  the  prob- 
lem. 

A.  simplified  model  of  electron  collision  frequency  as  a function 
of  electric  field  strength  is  given  which  enables  the  model  to  predict 
accurate  results  for  nuclear  weapon  gamma  yields  up  to  at  least  100  Kt. 
The  results  predicted  by  this  EMP  model  compare  to  within  5.5*  with 
results  from  the  Air  Force  Weapons  Laboratory  CHEMP  computer  code. 

The  computation  time  using  the  presented  model  on  a CDC  &600  compu- 
ter Is  typically  5 sec  or  less  for  a 5 shake  computation  period  in 
steps  of  .1  shake. 

The  mode!  presented  should  be  useful  for  both  classroom  instruc- 
tion and  nuclear  vulnerability/survivability  studies  and  analysis  prob- 
lems. 
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A CALCULATIOHAL  MODEL  FOR  HIGH  ALTITUDE  EHP 
I . Introduction 

In  recent  years  there  has  been  a growing  concern  about  the  nuclear 
weapon  effect  known  as  the  electromaonetlc  pulse  or  the  EMP.  As  the 
theory  of  EMP  developed  It  was  realized  that  there  are  actually  several 
types  of  EMP  which  are  characterized  by  the  mechanisms  which  produce 
them.  Types  of  EMP  Include  surface  EMP,  system  generated  EMP  (SGEMP) , 
and  high  altitude  EMP.  A comprehensive  discussion  of  the  va  'ous  types 
of  EHP  Is  presented  by  Kinsley  (Ref  I).  This  report  will  address  only 
one  of  these,  that  of  high  altitude  EHP. 

The  EMP  Is  basically  a long  range  nuclear  weapons  effect  wherein 
a high  Intensity  electromagnetic  field  Is  radiated  over  a wide  frequency 
band.  The  duration  of  this  pulse  Is  typically  on  the  order  of  shakes 
(1  shake  » I0~8  sec).  Specifically,  high  altitude  EMP  Is  that  produced 
by  a nuclear  detonation  above  20  Km. 

The  high  altitude  EMP  problem  is  a subject  of  great  Interest  to 
the  USAF  due  to  its  long  range  nature,  and  much  effort  goes  Into  pre- 
dicting the  time  development  of  the  fields.  These  high  Intensity  EMP 
fields  present  a potential  threat  to  weapons  systems  which  rely  on 
electrical  or  electronic  components,  such  as  communications,  electronic 
counter  measures  (ECM),  navigation,  guidance,  reconnaissance,  and  many 
others.  The  threat  levels  predicted  by  theoretical  models  have  a di- 
rect Impact  on  the  design  of  future  Air  Force  weapons  systems.  These 
threat  levels  also  help  to  predict  the  survivability  of  present  weapons 
systems. 
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Th«>  EMP  generation  models  currently  In  use  by  the  Air  Force  Wea- 
pons Laboratory  (AFWL),  are  highly  sophisticated  and  use  large  amounts 
of  computer  time.  Other  models  available  are  very  simple  and  give  es- 
sentially order  of  magnitude  results.  The  mode'  presented  in  this  re- 
port Is  intended  to  provide  an  alternative  betwetn  these  extremes, 
giving  results  close  to  those  of  the  more  advanced  AFWL  models  but  with 
very  short  expenditures  of  computer  time.  One  important  capability  of 
this  model  is  its  ability  to  account  for  a preionization  level  due  to  a 
precursor  burst. 

This  model  should  be  useful  for  quickly  obtaining  a meaningful 
estimate  of  the  EMP  environment  for  parameter  variation  or  for  use  in  a 
classroom  situation  to  give  a feel  for  the  calculations  involved  In  pre- 
dicting the  generated  EMP  fields  and  showing  the  effects  of  a given 
parameter. 

The  theory  behind  the  high  altitude  EMP  model  presented  is  based 
on  that  of  Karzas  and  Latter  (Ref  2).  Basically  the  Karzas-latter 
theory  states  that  prompt  gamma  rays  from  the  weapon  produce  Compton 
electrons  within  a specified  region  of  the  atmosphere  known  as  the  ab- 
sorption region.  These  Compton  electrons  are  turned  by  the  magnetic 
field  of  the  earth  to  produce  the  radiated  electromagnetic  fields. 

There  are  several  assumptions  which  are  basic  to  the  Karzas-latter 
model.  These  assumptions  are:  the  earth's  magnetic  field  is  assumed 

to  be  uniform,  the  magnetic  field  lines  have  no  curvature  in  the  gamma 
ray  absorption  layer,  the  earth's  surface  Is  assumed  to  be  flat,  the 
gamma  rays  are  monoenerglc,  each  gamma  Interacts  with  the  atmosphere  to 
produce  only  one  forward  directed  Compton  electron,  the  Compton  elec- 
trons produced  are  also  considered  to  be  monoenergic  and  have  a constant 


2 


gep/PH/75-13 

velocity  throughout  their  lifetime,  and  only  the  high  frequency  portion 
of  the  pulse  is  considered. 

One  additional  approximation  Is  made  In  this  study.  This  is  a 
small  angle  approximation  in  the  trigonometric  expressions  for  the  Comp- 
ton currents  and  air  conductivity..  This  approximation,  with  the  judi- 
cious choice  of  a function  to  represent  the  time  dependence  of  the  wea- 
pon yield,  leads  to  closed  form  expressions  for  the  Compton  currents 
and  air  conductivity.  This  saves  considerable  computation  time  and 
greatly  simplifies  the  calculations.  For  this  approximation  to  be 
valid  the  cyclotron  frequency  u of  the  Compton  electrons  In  the  qeonag- 
netic  field  must  be  small  enough  so  that  the  small  angle  approximations 
for  sin  ot  and  cos  wt,  where  t Is  in  seconds,  are  valid. 

In  all  cases  where  quantities  vary  with  the  atmospheric  density, 
an  exponential  atmosphere  model  is  used. 

Factors  not  considered  In  the  model  are  recombination  of  ions, 
avalanching  or  cascading  of  secondary  electrons.  X-ray  effects,  and 
self-consistent  electromagnetic  fields. 
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1 1 . Theory 


Overview 

Since  the  solution  to  the  EMP  problem  is  actually  the  solution  of 
a classical  electromagnetic  theory  problem,  the  derivation  of  the  model 
equations  reduces  to  putting  Maxwell's  equations  into  a convenient  form. 
This  is  essentially  accomplished  by  expressing  Maxwell's  equations  in 
spherical  coordinates  and  transforming  to  a retarded  time  frame.  One 
must  also  develop  expressions  for  the  currents  and  conductivities  of 
the  system  in  the  absorption  region.  The  general  equations  describing 
the  high  altitude  model  of  Karzas  and  Latter  have  been  derived  in  great 
detail  by  Chapman  (Ref  3).  Only  the  major  points  of  the  derivation  will 
be  given  here.  The  system  origin  is  assumed  to  be  at  the  burst  point 
with  detonation  at  time  t = 0.  This  geometry  is  illustrated  in  Fig.  I. 

The  key  points  to  be  remembered  in  this  model  are: 

1.  Each  gamma  ray  gives  rise  to  one  downward  traveling 
Compton  electron. 

2.  The  electrons  are  turned  by  the  earth's  magnetic  field 
giving  rise  to  a centrifugal  acceleration. 

3.  The  relativistic  electrons  radiate  energy  in  their 
forward  direction. 

A.  The  gamma  rays  and  the  EMP  radiation  travel  at  the 
same  speed.  This  leads  to  constructive  interference 
of  the  radiation  from  each  of  the  electrons. 

Particle  Densities 

The  gamma  rays  from  >.ne  nuclear  weapon  travel  in  a straight  line 
to  a point  where  they  produce  Compton  electrons.  At  any  given  point  r 
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t’ic  number  of  gammns  which  interact  to  produce  Compton  electrons  Is 


expf-/,4C] 
v o *(r  > 

£"  *orr2  X (r) 


whore  X(r)  is  the  mean  free  path  of  gamma  rays  to  produce  Compton  elec- 
trons, Y is  the  gamma  yield  of  the  weapon  ir.  electron  volts  (eV),  and  E 
is  the  mean  gamma  energy  in  eV, 

Equation  (l)  may  also  be  called  the  radi3l  distribution  function, 

Y 

or  an  attenuation  function  for  interacting  gamma  rays.  The  g-  term  is 
the  total  number  of  gamma  rays  available  from  the  weapon.  The  4«r2 
term  accounts  for  the  divergence  of  the  gamma  rays  as  the  radius  r Is 
Increased  while  the  remaining  terms  account  for  the  reduction  In  gammas 
due  to  th-.ir  absorption  in  the  atmosphere,  based  on  the  mean  free  path. 

It  is  assumed  that  the  gamma  mean  free  path  varies  as  the  exponen- 
tial atmosphere.  This  gives  the  functional  relationship  between  \ and 


Mr)  * X g exp  [(HOB  - r cos  A),-„ 


X q « gamma  mean  free  path  at  standard  pressure 
HOB  « height  of  burst  in  Km  above  the  earth's  surface 
r ■ radial  distance  from  the  burst  point  to  the  point  of  Interest 
A « angle  between  the  position  vector  r and  the  vertical 
S » atmospheric  scale  height 

With  this  assumption,  Eq  (1)  can  be  Integrated  and  becomes 


, , Y 1 f S 

9 (r)  - r T — TTTT  exp  i * , 7 

c liirr2  X(r)  **  Aq  cos  A 


r HOBn  r f r cos  A*\  ,1 1 

exp  —J  |_exp  — $ — ) - 1 J J 
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Mow  if  f(t)  is  the  time  distribution  function  of  the  weapon  yield, 
the  rate  of  Compton  electrons,  nc,  produced  at  a given  point  r and  time 


t is  given  by 


r-<«  '0-0 


Each  Compton  electron  produces  through  inelastic  scattering  events 
several  secondary  electrons  which  form  the  basis  for  the  conductivity 
of  the  atmos- '•ere.  As  in  the  Kar2as-Latter  approach,  each  Compton  elec- 
tron is  assumed  to  have  a constant  speed,  VD,  throughout  the  ranoe,  R, 
of  the  electron  which  is  a function  of  altitude.  This  allows  the  life- 
time to  be  expressed  as  R/V0.  If  each  Compton  electron  produces  secon- 
dary electrons  at  a constant  rate,  the  rate  of  secondary  electron,  ns, 


production  is 


dns  Ec/33  «V 


where  Ec  is  the  energy  of  the  Compton  electron  and  33  eV  is  the  average 
ionization  energy  per  air  molecule  (Ref  2). 

Considering  the  differential  current  produced  by  the  Compton  elec- 
trons, it  can  be  shown  (Ref  3)  that  the  Compton  current  and  the  number 
of  Compton  electrons  are  given  by 


Jc(t)  = - e g(r)  J ° V(t')  f (x  - t'  + * — -*)  dt' 

O C y 


nc(t)  = g (r)  / ° f (r  - t'  + dx' 


-v  ^ 


fit P/PH/7 5- 13 


i'  “ the  time  since  the  creation  of  the  Compton  electron 
X(t")  ■ the  radial  distance  the  Compton  electron  has  traveled 
e « the  magnitude  of  the  electron  charge 
The  quantity  t is  generally  known  as  retarded  time. 

It  then  follows  from  Eqs  (A)  and  (7)  that  the  number  of  secondary 
electrons  is 


r r°  f G'  - f ^ 


X(t")>  . ..1  . . 

J dt  ' I dt 


where  q is  Ec/33  eV. 


Currents  and  Conductivity 

In  the  Karzas-Latter  theory,  the  speed  of  the  Co-ipton  electrons  Is 
considered  to  be  a constant,  however,  there  is  an  acceleration  due  to 
the  geomagnetic  field.  The  general  equation  of  notion  for  an  electron 
in  this  case  is 


L.  n yV  = - e (E  + V X B)  - ra  vcV 
dt 


"»  = the  electron  rest  mass 

V = the  electron  velocity 
E * the  electric  field 

B * the  nagr.etlc  field 
vc  * the  electron  collision  frequency 

Y “ (1  * (V0/c)2)*^ 
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If  the  relativistic  motion  of  the  electron  is  considered,  only  the 
V X B terra  is  iranortont  and  in  sohcric.il  coordinates,  the  expressions 
for  the  velocity  components  become  (Ref  3) 


Vr  = V0  (sin2  6 cos  ot  + cos2  0) 


V,  -■  V„  (cos  9 sin  9 cos  ot  - sin  9 cos  9)  (11) 

O O 


Vg  = V0  (sin  9 sin  ot) 

where  o i"  the  cyclotron  frequency  for  an  electron  and  is  given  by 


with  0o  the  magnitude  of  the  geomagnetic  field. 
From  Eq  (10),  X(t')  is  found  to  be 


X(t')  = V0  f sin2  9 S '-nj-T—  + t'  cos2  bA 


The  Compton  currents  may  nov;  be  written  as 


Jc(t)  - - cg(r)  V0  / ° (f(T)(cos2  9 + sin2  9 cos  ot)l  dt'  (15) 


R/V. 

JC(t)  * - eg (r)  V / [f  (T)  sin  6 cos  8 (cos  ut'  - D]  dt' 


jc(t)  » - eg(r)  V / ° tf(T)  sin  9 sin  ot']  dt' 
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where  ft  = — . 
c 

In  a like  manner  £q  (20)  becomes 

T'  - t'  - (1  - 6)  t"  (77) 


An  expression  for  the  conductivity  may  also  be  found  by  using  the 

equation  of  motion  for  the  secondary  electrons.  These  electrons  are  in 

the  thermal  regions  with  energies  ranging  from  about  10-15  eV  to  the 

amoienc  energy.  It  should  be  remembered  for  later  use  that  the  ambient 

energy  of  the  secondary  electrons  is  dependent  on  the  electric  field 

present.  For  consideration  here,  it  is  assumed  that  y \ I and  also  that 

the  change  of  velocity  with  time  is  small  compared  to  the  other  terms 

in  Ea  (9)  so  that  nay  be  neglected.  Also,  with  low  velocities,  the 

dt 

V X I term  is  small  compared  to  the  remaining  terms  and  may  also  be 
neglected.  Then  the  velocity  of  the  secondary  electrons  is 

V „ - S_  E (28) 

mvc 

Using  Eq  (28),  the  current  due  to  the  secondary  electrons  is 

Js(x)  - - e V ns(T>  = ns(t)  (29) 

Comparing  Eq  (29)  to  Js  » oE,  an  expression  for  the  conductivity 
is 

0(t).sL„s(t)  (30) 

•'m'e 

Equations  (19),  (23).  (2*t),  (25),  (26).  (27),  and  (30)  provide  the 
desired  expressions  for  the  Compton  currents  and  the  conductivity. 
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Field  Equot Ions 

Maxwell's  equations  is  rational i zed  MKS  units  are 


•*  •+  aB 

V X E ® • ~ 
at 


•*  + | ££ 
7 X B = ucJ  + 


■*■  * Qv 

7.£  - — 

£0 


where  qv  is  the  total  charge  density  and  J is  the  total  current  density. 
In  addition  to  these  equations  the  continuity  of  charge  requires  that 


Combining  these  equations  to  separate  E and  E and  transforming 

them  into  spherical  coordinates  and  into  the  retarded  tine  frame  (Ref 
-♦  •* 

3)  the  relations  for  E and  S become 


i ♦ ♦ l -*■ 

72E  + ur  - — V.J  + — Vq„ 
r «o  c0  V 

3 7 2 1 3 , *,  - . . n 

+ 57  .77  57  (rE)  + uo  (J  * ur  Jr>  j 


W-Po’*J*57  [ 77  57  (rB) 


+ 7-  - U0V  J “ 0 


In  the  Karzas-Latter  model,  only  the  time  derivative  portion  of 
Eqs  (36)  and  (37)  are  kept  since  the  current  variation  with  distance 
is  slow  compared  to  the  variation  in  tine  for  the  high  frequency  com- 
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ponents.  Also  the  fields  and  currents  vary  r op  idly  in  time.  This 
approximation  is  valid  for  about  100  shakes.  Th  s sane  hioh  f reoucncy 
approx inat ion  is  used  here.  In  addition,  the  radial  component  of  the 
field  is  dropped  since  it  is  weak  compared  to  the  transverse  components 
and  contributes  only  a very  low  frequency  sionol  (Ref  2).  The  equa- 
tions for  the  transverse  comnonents  are 


It  T f fIf  (rEe.o)  + uo  J«.f] 

(38) 

(39) 

k'i'rJrW  *?*]•' 

(90) 

The  currents  in  Eqs  (33),  (30),  and  (^0)  are  total  currents.  The 
total  currents  are  given  by 

Je,e  “ J0,2  + o(t)  E8,0 

(91) 

Substitution  of  Eq  (91)  into  Eqs  (33),  (39),  and  (90)  and 
tion  over  time  gives 

intepra- 

I 7 57  (rEj  + u5Jg  + y0  c(t)  Eg  = 0 

Iff?  (rV  + BOJ0  + **0  3(t>  EfJ  - 0 


I?! Es>'° 


Iff?  (pty  + -J9  + ~o(t)  E6  ' 0 


m 

m 

m 


('.5) 
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Equations  (f»?)  and  (^3)  are  In  a form  which  can  be  solved.  The 
ternt  needed  for  solution  of  these  equations  arc  the  air  conductivity 
and  the  transverse  components  of  the  Compton  currents.  The  air  con- 
ductivity nay  be  found  by  using  Eqs  (19)  and  (30)*  The  transverse  com- 
oonents  of  the  Conpton  currents  may  be  found  by  using  Eqs  (2*»)  and  (25). 
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ill.  MoJcl  Cons Herat ions 


While  the  basic  relationship  expressions,  have  been  developed  be- 
tween  E,  J,  and  o,  sene  additional  considerations  are  necessary  in 
develooino  the  presented  model.  These  include  such  items  as  the  form 
of  the  function  f(t),  the  average  energies  of  the  qamna  rays  and  Comp- 
ton electrons,  the  mean  free  path  of  the  na.mma  rays  and  the  range  of 
the  Compton  electrons,  and  the  collision  frequency  of  the  secondary 
electrons.  The  validity  of  the  small  angle  approximation  must  also  be 
examined.  Each  of  these  will  be  discussed  separately  in  this  section. 

When  all  of  the  items  discussed  in  this  section  a-?  combined  with 
the  equations  given  in  the  previous  section,  one  obtains  a straight- 
forward ard  simple  calculat ional  model  for  predicting  EMP  environments 
from  high  altitude  nuclear  bursts. 


Time  Dependence  of  Gamma  Rays 

Many  functions  have  been  proposed  to  describe  the  time  dependent 
output  of  a nuclear  weapon,  f ( t ) • These  pulse  shapes  are  often  taken 
to  be  of  the  form  (Ref  *») 


f(t)  = exp  (-  ?-  - St)  U(t) 


where  U(t)  is  a unit  step  function  havina  a value  zero  for  negative 
arguments  and  unity  for  ail  others,  and  a and  6 are  constants  describ- 
ing the  rise  and  fall  of  the  pulse.  Another  form  often  used  (Ref  5)  Is 


(q  + S)  exp  n (t  - tp) 

6 + a exp  [ (a  + S)  (t  - t0)] 
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where  t0  is  the  position  of  the  peak  value  of  the  function.  One  fea- 
ture of  these  functions  is  that  they  cannot  be  intenrated  :n  closed 
form,  and  thus  EMP  codes  which  use  the  functions  involve  ime  consuminq 
numerical  integration  schemes  for  calculating  the  Compt  currents  and 
the  atmospheric  conductivity. 

A basic  idea  in  this  investigation  was  to  select  functional  forms 
of  f(t)  that  could  be  integrated  directly,  thus  yielding  substantial 
savings  in  computation  time.  There  are  two  constraints  that  the  func- 
tional form  of  f(t)  should  follow.  The  first  constraint  is  that  the 
function  be  integrable  in  closed  form.  A second  constraint  is  that  the 
function  should  match  reasonably  well  those  currently  in  use  for  EMP 
calculations. 

Since  each  of  the  functions  given  in  Eos  {k&)  and  (A7)  has  an 
exponential  rise  and  decay,  prime  consideration  was  given  to  functions 
having  these  same  characteristics.  With  this  in  mind,  the  field  was 
narrowed  down  to  the  following  normalized  functions: 

f(t)  » lexp  (-<*0  - exp  (-St)]  U(t)  (A8) 

& - o 

...  l [exp  (St)  - 1]  U(t),  t <.  t 

f(t)=N 

[exp  <Bt0)  - I]  [exp  - a(t  - t„)]  U(t  - t0)  (1>9) 

where  H is  a normalization  constant. 

It  should  be  apparent  at  this  point  that  the  first  constraint 
given  could  not  be  met  without  the  simplification  made  to  Eqs  (16), 
(17).  08),  and  (20)  by  the  small  angle  approximations  In  Eqs  (21)  and 
(22). 
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After  fui  ther  study,  it  was  determined  that  Cq  (43)  would  fit  the 
functions  qiven  in  Eqs  (46)  and  (47)  closely  except  for  very  narrow 

pulses.  Figure  2 shews  the  shape  of  a typical  long  pulse.  It  is  a plot 

of  Cq  (4£)  with  a = I07  and  6 <=  3.7  (10)8.  The  X narks  indicate  the 
values  computed  from  Eq  (46)  with  a = I0"8  and  6 » 107.  Figure  3 shows 

a plot  of  a typical  narrow  pulse.  In  this  case,  Eq  (A3)  with  a = 2 

(10)8  and  S = 3 (10)8,  and  Eq  (46)  with  a « 2 (10)*8  and  8 ■ 2 (I0)8 

viere  used.  The  values  fron  Eq  (46)  are  again  annotated  by  X marks.  The 

curve  computed  from  Eq  (48)  has  been  shifted  in  time  away  fron  the  ori- 
gin in  order  for  the  peaks  of  the  two  functions  to  match.  This  was  done 

for  an  easier  comparison.  The  result  of  this  manipulation  is  to  offset 

the  times  in  the  final  calculations  by  the  amount  necessary  to  shift 
Eq  (48).  In  both  Figs.  2 and  3,  the  functions  were  normalised  to  a 
peak  value  of  one  for  comparison  purposes. 

Equivalent  comparisons  using  Eq  (49)  in  place  of  Eq  (48)  indicated 
that  Eq  (43)  would  give  a better  overall  versatility.  Thus  Eq  (48)  was 
chosen  as  the  pulse  shape  fur  tne  sample  calculations  of  this  E.’IP  model. 


Energies 

Two  energy  values  are  needed  in  the  model,  that  of  the  gamma  rays 
and  that  of  the  Compton  electrons.  Karzas  and  Latter  (kef  2)  suggest 
a I t'.cV  gamma  and  a ! MeV  Compton  recoil  energy.  This  implies  no 
energy  loss  in  the  Compton  creation  process.  Chapman  (kef  3)  uses  1.5 
MeV  gamma  rays  and  I.7.S  HeV  Compton  electrons,  which  are  the  most  ener- 
getic available  from  1.5  HeV  gamma  rays.  The  values  used  in  the  AFWL 
models  vary  somewhat  (Ref  4),  but  the  most  common  value  appears  to  be 
1.5  HeV  gamma  rays  and  .75  HeV  Compton  electrons.  These  Comoton  elec- 
tron energies  are  actually  the  kinetic  energy  transferred  to  the  elec- 
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tron  from  the  gamma  ray  and  assume  the  Compton  electrons  are  only  scat- 
tered forward.  The  encroy  of  .75  MeV  from  a 1.5  HeV  qanma  ray  corre- 
sponds to  the  averanc  encroy  transfer  in  a Compton  creation  process 
(Ref  6).  This  Compton  electron  enerqy  is  qiven  by 

T = (1  - as/oc)  hv  (50) 

where 

T = average  recoil  energy 

os  «■  average  Compton  scatterinq  (non-absorption)  cross  section 
oc  = average  Compton  collision  (total  interaction)  cross  section 
hv  = energy  of  the  gamma  ray 

The  values  of  as  and  ac  are  from  the  Klein-liishina  formulas  (Ref  6). 

With  these  considerations,  1.5  HeV  gamma  rays  are  assumed  while 
the  Compton  energy  is  left  as  a variable. 

A question  directly  related  to  the  energies  is  that  of  gamma  mean 
free  path  and  Compton  range.  The  mean  free  path  is  defined  as  the  in- 
verse of  the  total  cross  section.  For  Compton  creation  the  total  macro- 
scopic cross  section  is  the  number  of  electrons  times  the  absorption 
microscopic  cross  section.  The  absorption  cross  section  is  the  differ- 
ence between  the  collision  (total  interaction)  and  scatterinq  cross 
sections.  Using  this  definition,  the  mean  free  path  varies  little  for 
gamma  energies  between  1 and  1.5  HeV  (Ref  6)  so  that  the  value  sugges- 
ted by  Karzas  and  Latter  of  300  n at  standard  pressure  Is  a realistic 
value  and  is  used  in  the  model. 

Katz  and  Penfold  (Ref  7)  give  relatlonshios  between  electron 
energy  and  range  in  any  material.  For  electrons  with  energies  less 
than  2.5  HeV  the  relationship  is 
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R = [412  E <!‘265  ‘ °-C95''  ,n  Eo>]/0 


(51) 


■ i 

i 


where  R is  the  ranqe  in  cm,  E0  is  the  energy  in  MeV  of  the  electron, 
and  0 is  the  density  in  mg/cm3  of  the  material  of  interest.  For  a .75 
MeV  electron  and  an  air  density  of  1.293  mq/cm3  at  standard  pressure, 
the  range  is  2.23  n.  For  a 1 MeV  electron  the  range  is  3-18  m,  which 
shows  an  increase  of  almost  one  meter.  Since  the  Compton  recoil  energy 
is  taken  as  a variable  and  the  range  variation  is  considerable,  Eq  (51) 
with  an  air  density  of  1.293  mg/cm3  is  used  in  the  model  to  compute  the 
range  of  the  Compton  electrons  at  standard  pressure. 

Another  factor  which  is  a function  of  Compton  electron  energy  is 
the  velocity  of  the  Compton  electrons.  The.  energy  used  for  the  Compton 
electron  is  considered  to  be  the  kinetic  energy  of  the  electron.  The 
total  energy  is  the  sum  of  the  kinetic  and  rest  energies.  The  relati- 
vistic relation  here  is 


1 


E + ro0  » m0y  (52) 

where  E is  the  kinetic  energy  in  eV  and  m0  is  the  rest  energy  in  eV. 
Equation  (52)  may  be  solved  for  the  Compton  electron  velocity  to  get 


Vo 


c 


mo 

E + mD 


(53) 


Equation  (53)  is  also  used  in  the  model. 

The  energy  from  the  prompt  garma  rays  is  assumed  to  be  deposited 
in  an  absorption  region  from  20  Km  to  50  Km  above  the  earth's  surface. 
This  is  based  on  calculations  by  Latter  and  lelevier  (Ref  8). 
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Compton  Electron  Lifetime 

The  lifetime  of  a Compton  electron,  as  previously  defined,  is 
given  by  the  quotient  of  the  electron  ranoe  P.  and  its  constant  speed 
V0.  The  ranoe  at  any  given  altitude  is  found  by  exoonent ia I ly  scaling 
the  standard  pressure  value  found  from  Eq  (51)  and  the  sea  level  air 
density.  Since  the  high  frequency  approximation  is  only  valid  for 
about  100  shakes,  the  lifetime  should  not  be  allowed  to  exceed  this 
length  cf  time. 

Small  Annie  Approximation 

One  inportant  factor  for  this  model  is  the  validity  of  the  small 
angle  approximations  leading  to  Eqs  (2f» ) , (25).  (26),  and  (27).  This 
can  be  somewhat  arbitrary,  depending  on  the  amount  of  error  one  Is 
willing  to  accept.  For  the  use  of  this  analysis,  a 1%  error  will  be 
considered  acceptable.  With  this  limit  the  maximum  value  for  the  solu- 

i 1 

tion  of  sin  A ■ A is  A = .36.  i i 

The  maximum  geomagnetic  field  in  the  qar.ma  ray  absorption  region  " 

is  usually  taken  to  be  a .6  gauss  (Ref  4).  A minimum  kinetic  energy  ; r 

J ; 

of  .5  MeV  will  be  assumed  for  the  Compton  electrons.  The  time  which  . . 

satisfies  : 's 

or  = .yh  (54) 

is  the  maximum  time  duration  over  which  the  small  angle  approximation 

is  val id.  ; 

Equation  (13)  gives  the  expression  for  j.  This  will  be  a ma  :imum 
when  B0  is  a maximum  and  y is  a minimum.  For  y to  be  a minimum,  the 
particie  speed  must  be  a minimum.  The  speed  of  a .5  MeV  electron  is  j 

I v 

.863c.  With  these  factors  u is  5.33  (I0)6.  Solving  for  t in  Eq  (5b),  | : 
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the  result  is  t ° 6.*i  shakes.  This  is  however  a very  worst  case.  In 
(rest  situations  the  electron  energy  will  be  considerably  greater  than 
•5  MeV.  For  on  electron  of  .75  McV  with  a speed  of  .9l**2c,  the  time 
limit  is  t s 7.9  shakes.  If  the  geomagnetic  field  Is  .3  gauss,  which 
is  often  the  case,  the  time  limits  above  will  be  doubled. 

The  minimum  length  of  the  confidence  interval  should  be  C.k 
shakes.  As  the  electron  energy  increases,  the  length  of  the  confidence 
interval  will  also  increase. 


Electron  Collision  Frequency 

In  previous  EH°  models  two  general  approaches  have  been  taken  to 
find  a value  for  the  air  conductivity.  The  first  uses  Eq  (30)  and 
assumes  a constant  collision  frequency  scaled  to  an  exponential  atmo- 
sphere. This  approach  tends  to  limit  the  usefulness  to  a rather  nar- 
row range  of  gamma  yields  if  accurate  answers  are  desired.  This  basic 
approach  is  the  one  suggested  by  Karzas  and  Latter. 

The  second  approach  is  the  use  of  the  electron  transport  equations 
(Ref  9)  to  solve  for  the  secondary  electron  velocity.  This  velocity  is 
then  used  to  compute  the  secondary  currents  directly.  These  equations 
however,  require  an  iterative  solution  which  is  very  time  consuming. 

Since  the  model  developed  here  is  desired  to  be  as  simple  as 
possible  to  use  and  understand,  but  still  give  good  results  over  as 
wide  a range  of  gamma  yields  as  possible,  an  intermed iate  approach  Is 
desirable. 

When  a secondary  electron  is  created,  it  has  an  energy  somewhat 
higher  than  Its  ambient  energy.  The  secondary  electrons  transfer  this 
additional  energy  to  air  molecules  through  both  elastic  and  inelastic 
collisions.  The  time  required  for  this  process  to  take  place  Is  known 
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as  the  thermal  izat ion  time.  Eaun  (Rof  10)  has  mothered  data  to  show 
that  toe  thermal  izat ion  time  for  low  energy  (orou?  10  cV  or  less)  elec- 
trons at  sea  level  is  approximately  0 nanosecond.  Us  inn  on  exponential 
atmosphere  model,  this  thermal i zat ion  time  corresponds  to  approximately 
2 shakes  at  20  Km  and  to  over  a microsecond  at  50  Km.  The  thermal  iza- 
tion  time  in  the  20  Km  case  is  on  the  same  order  as  the  time  scale  of 
interest  (a  few  shakes)  for  the  high  altitude  EM*5  problem.  The  scale 
height  used  for  the  exponential  atmosphere  is  taken  as  7 Km.  This  value 
is  based  on  a curve  fit  of  data  from  the  U.S.  Standard  Atmosphere  of 
1 962  (Rer  II).  This  scale  height  of  7 Km  is  used  throuohout  the  model. 

Baum  also  shows  the  effect  of  an  applied  electric  field  on  the 
ambient  electron  energy.  The  secondary  electron  ambient  energy  in- 
creases with  increasing  field  strengths  above  10H  volts/m.  For  elec- 
tric field  strengths  below  104  volts/m  there  is  little  or  no  effect. 

If  the  thermal  ization  oroccss  is  assumed  to  be  a linear  function 
of  tine  and  the  lower  energy  limit  of  the  secondary  electrons  is  as- 
sumed to  be  a linear  function  of  the  electric  field  strength,  these 
effects  can  be  easily  incorporated  into  the  EMP  rode  I be  inn  developed. 

It  is  however,  necessary  to  determine  the  coefficients  of  the  assumed 
functions. 

The  tine  dependent  function  is  based  on  the  thermal izat ion  tine 
at  20  Kn.  The  coefficients  for  the  electric  field  dependent  function 
are  based  on  an  empirical  fit  to  the  peak  electric  field  values  from 
the  AFWL  CHEMP  model  (Ref  A).  In  order  to  get  an  acccntable  fit,  It 
was  found  that  the  field  dependent  function  had  to  be  split  into  two 
linear  functions.  Since  the  electron  energy  and  collision  frequency 
are  directly  related  (Ref  9),  the  linear  relations  described  above  may 
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be  expressed  directly  as  electron  collision  frequency.  The  functions 
described  arc 


-2.5  do)20  T + k.5  (I0)12 

(55a) 

0.k3  (10)s  E + 1.6  (I0)12 

E < s do)1* 

(55b) 

6 (10)7  E + 0.6  (I0)12 

E > 5 (10)" 

(55c) 

where  t is  retarded  tine  in  seconds  and  E is  the  electric  field 
strength  in  volts  per  neter. 

The  collision  freouency  at  sea  level  is  taken  as  the  maximum  of 
Eos  (55)  and  2.S  ( 1 0) 1 2 . In  no  case  is  the  collision  frequency  allowed 

to  exceed  k.k  {1 0) 1 2 . The  collision  freouency  is  then  scaled  using  an 
exponential  atmosohere.  The  lower  limit  value  of  2.3  (10)12  was  deter- 
mined empirically  while  the  upocr  limit  value  of  k.k  (lo)12  results 
from  the  collision  frequency  being  essentially  independent  of  the  elec- 
tric field  for  extremely  high  fields  (Ref  2). 

Preion i zat ion 

The  problem  of  preionization  by  a precursor  burst  may  be  handled 
quite  easily  by  the  model  developed  here.  This  is  because  of  the 
choice  of  the  oamna  yield  function  of  Eq  (kC)  and  the  small  angle  ap- 
proximations of  Eqs  (21)  and  (22),  which  allows  Eq  (I?)  to  be  Inteqra- 
ted  in  closed  form  to  get  the  nunber  of  secondary  electrons  as  a func- 
tion of  retarded  time  for  any  altitude.  The  total  air  conductivity 
for  a burst  with  a precursor  burst  is  comouted  from  the  total  number 
of  secondaries  from  both  bursts  at  the  time  of  interest.  The  colli- 
sion frequency  should  be  based  on  the  time  since  the  precursor  burst 
and  the  electric  field  strength  due  to  the  main  burst  at  the  time  of 
interest. 
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it  is  In  the  calculation  of  the  number  of  secondary  electrons  from 

a precursor  burst  where  the  tine  confidence  limit  of  the  snail  angle 

approximations  is  greatly  exceeded.  The  term  of  interest  here  is 

LI — , — in  Eq  (20).  Since  this  term  and  its  small  amle  approximation 

u 

t"  both  have  maximum  on  the  order  of  I0"6,  and  recombination  Is  a 
noticeable  factor,  on  preionization  time  scales,  any  error  introduced 
by  the  small  angle  approximation  should  be  tolerable. 


Calculation  Method 

The  time  development  of  the  EKP  signal  may  be  found  by  the  numeri- 
cal integration  of  Eqs  (A2)  and  (A3)  over  the  absorption  region  along 
the  line  of  sight  between  the  burst  point  and  the  target  for  retarded 
times  from  0 to  i.  The  distance  from  the  burst  point  to  the  top  of 
the  absorption  region  along  the  line  of  sight  is  called  RMIti  while  the 
corresponding  distance  to  the  bottom  of  the  absorption  region  is  called 
RMAX.  The  distance  to  the  target  from  the  burst  point  is  r.  It  should 
be  noted  that  above  and  below  the  absorption  region  there  are  no  cur- 
rents or  air  conductivity  so  that  the  only  contribution  to  the  electric 
field  is  over  the  absorption  region.  The  integration  limits  are  then 
RM1N  and  RHAX  as  shown  in  Fig.  I on  page  5.  RMAX  will  equal  r If  the 
target  is  in  the  absorption  region. 

Because  a numerical  integration  Is  necessary,  this  model  is  best 
run  on  a computer.  The  results  given  in  this  report  were  calculated 
through  a computer  code  usino  the  concepts  and  methods  detailed  here 
applied  to  the  theory  given  previously.  The  computer  code  calculation 
will  be  referred  to  as  HAEMP  in  the  next  section. 

For  points  below  the  absorption  region,  where  both  the  Compton 


currents  and  air  conductivity  are  zero,  Eqs  (A2)  and  (A3)  have  soiu- 
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hold,  the  first  term  in  Eqs  (42)  and  (43)  may  be  neglected  and  they 


Eg  » - Jg/o(r)  (61 

Eg»-J§/o(t)  (62 

the  values  found  from  Eqs  (6l)  and  (62)  may  then  be  used  with  Eq  (53) 
to  find  the  electric  field  at  the  target. 
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IV,  Pesul ts 

To  test  the  validity  of  the  liAEMP  model,  the  ccmnuted  values  of 
this  model  were  compared  to  available  equivalent  model  values  from  the 
AFWL  code  for  EMP  calculations  t.nowr  as  CIIEttP  (P.ef  A). 

The  basic  set  of  conditions  used  for  these  calculations  was: 
target  location  = nround  zero 

height  of  burst  = 100  Km 

geomagnetic  field  = .3  nauss  or  3 ( « 0) ~s  wb/m2 

inclination  angle  = 0.0  degree 

Comptcn  electron  recoil  energy  = .75  HeV 

The  gamma  yield  was  va  :ed  and  the  number  of  steps  taken  for  the 
numerical  integration  was  varied  according  to  the  gamma  yield,  with 
more  steps  taken  for  the  higher  yields.  Other  parameters  varied  for 
examination  of  ocak  field  values  were  burst  height,  geomagnetic  field, 
and  pulse  shape.  A preioni zat ion  level  was  also  considered. 

The  available  data  from  the  CHEMP(ll)  code,  which  is  CHEMP  run  with 
non-self-consistent  calculations,  was  computed  using  a pulse  of  the 
form  of  Eg  (46)  with  a » 10*8  and  S = 1C7  and  the  same  geometry  as 
given  above.  This  pulse  shaoe  is  almost  identical  to  that  of  Eq  (AS) 
with  a 107  and  6 = 3.7  (10)8.  See  Fig.  2 on  page  !£.  Using  this 
pulse  shape,  a range  of  gamma  yields  from  .01  Kt  to  100  Kt  was  used  to 
calculate  the  EMP  field  values.  The  peak  field  values  are  plotted  in 
Fiq.  4.  The  values  taken  from  the  CHEMP(N)  code  are  annotated  by  X 
marks.  These  peak  field  values  shew  a maximum  difference  of  5.5% 
around  the  I Kt  case,  and  a difference  of  less  than  2%  for  all  other 
known  cases. 
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The  complete  pulse  out  to  i < ihak.es  for  the  .2D  Kt  case  is  shown  In 
rjq.  5.  '■.gain,  values  from  the  '.HW  ('0  code  are  indicated  by  X marks. 
There  is  excellent  agreement  for  the  first  two  shakes  after  wh'ch  time 
the  IIAEKP  values  fall  off  more  rapidly.  This  more  rapid  fall  is  most 
likely  caused  by  two  factors.  The  first  is  not  considering  recombina- 
tion of  the  secondary  electrons,  which  would  tend  to  reduce  the  con- 
ductivity and  allow  the  fields  to  faP  off  more  slowly.  The  second  is 
the  slew  breakdown  of  the  small  angle  approximation  in  computing  the 
Compton  currents. 

The  height  of  burst  was  varied  to  observe  the  effect  on  the  peak 
fields  for  the  various  yields.  This  data  is  plotted  in  fig.  6 for 
altitudes  of  60  Krt,  100  Km,  200  Km,  and  300  Kn.  The  neneral  trend  of 
Fig.  6 indicates  that  the  field  strength  increases  as  the  height  of 
burst  increases  only  for  the  larger  yields,  for  the  lower  yields,  the 
field  strength  decreases  as  the  burst  height  increases.  Examination 
of  the  hypothetical  curves  in  fig.  6 implies  that  the  maximjm  EMP  sig- 
nal generated  by  any  given  yield  depends  on  the  height  of  burst,  for 
example,  using  the  curves  in  Fig.  6,  it  could  be  concluded  that  the 
maximum  EHP  from  a I Kt  gamma  yield  would  occur  for  a height  of  burs 
somewhere  between  60  Km  and  200  Kn  and  for  a . 1 Kt  oamma  yield,  the 
burst  height  would  be  less  than  100  Km. 

A test  was  also  run  to  determine  the  effect  of  a narrower  gemma 
output  pulse.  A pulse  used  by  the  AFVfl  in  the  CHEMP  code  for  this 
ourpose  is  of  the  form  of  Eq  (A6)  with  a - Z (10)  8 and  3 * 2(10)8. 

An  essentially  identical  pulse  can  be  achieved  with  Eq  (W)  using 
a * 2 (10)9  and  3 - 3 (10)8 . See  Fig.  3 on  page  19-  A plot  of  peak 
fields  vs  yield  Is  shown  in  Fig.  7 for  both  the  narrow  pulse  and  the 
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Fig.  6.  Peak  Electric  Field  as  a Function  of  Curst  Height 
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original  long  pulse.  This  data  points  to  n saturated  peak  field 
strcnath  of  b.32  (10)1*  volts/m  for  the  100  Yn  burst  case  in  a .3  oauss 
geomagnetic  field.  The  narrow  pulse  curve  shows  that  the  peak  fields 
ore  higher  for  all  yields.  The  data  from  Fin.  6 however  shows  that 
the  saturated  field  value  preoictcd  in  Fin.  7 does  not  hold  for  all 
burst  heights. 

Another  factor  which  should  have  a noticeable  effect  on  the  EMP 
generated  is  the  presence  of  a precursor  burst.  For  this  test  a .03 
Kt  qamna  yield  burst  was  assumed  to  have  bpen  detonated  10  nsec  prior 
to  the  main  burst.  The  pulse  shape  used  to  describe  the  qamma  output 
is  the  initial  long  pulse  ana  Is  used  for  both  bursts.  The  peak  field 
results  are  plotted  ir,  Fie,  G.  The  results  of  this  test  were  au.te 
dramatic.  This  data  shows  that  the  pear,  field  of  1.6  (ID)4  volts/n 
from  the  .03  Kt  precursor  burst  is  not  equaled  until  the  main  burst  has 
a yield  of  nearly  1 Kt.  In  general,  the  neak  field  values  with  this 
level  of  oreicni7ation  are  tho  same  as  If  the  weapon  yield  was  one 
tenth  of  its  actual  yield,  As  the  yield  of  the  main  burst  goes  up,  the 
field  reduction  becomes  less.  At  a yield  of  .1  Kt  the  field  is  reduced 
by  over  a factor  of  A.  At  .25  Kt  the  factor  drops  to  3 and  to  2 for 
the  I Kt  case.  When  the  yield  finally  reaches  100  Kt,  the  reduction 
factor  has  dropped  to  1.1. 

The  results  with  preionization  compare  favorably  with  available 
results  from  the  AFWL  (Ref  A).  However,  since  a slightly  different 
geometry  and  some  additional  source  modeling  were  used,  a direct  com- 
parison cannot  be  made. 

A last  major  factor  to  consider  for  this  model  is  the  geomagnetic 
field  strength.  The  peak  field  values  as  a function  of  yield  are 
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plotted  in  Fig.  3 for  the  .3  gauss  ana  the  .6  gauss  ocenannctic  fields. 
The  larger  geomagnetic  field  resulted  in  apnroxinately  a tines  2 In- 
crease in  the  field  strengths  for  all  vields  considered. 

One  renaining  result  should  be  mentioned.  This  concerns  the  tine 
of  the  peak  field  value.  Since  for  these  calculations  the  fields  were 
computed  in  steps  of  .1  shake,  the  peal  field  tine  will  he  accurate 
only  to  this  Unit.  The  peak  for  the  .25  Kt  case  shown  in  Fig.  5 oc- 
curs at  1.7  shakes  for  the  IIAEHP  node!  and  at  an  estinated  1.3  shakes 
for  the  CHEMP(N)  code.  As  the  yield  goes  up,  the  peak  field  value 
occurs  earlier.  For  yields  above  .25  Kt,  the  peak  occurs  between  1.1 
and  1.3  shakes  for  all  cases  examined.  For  all  yields  and  altitudes 
considered,  the  peak  always  occurred  within  8 shakes,  with  the  majority 
falling  between  1 and  2 shakes. 

Since  the  model  presented  is  supposed  to  be  a quick  computational 
tool,  some  estimates  on  the  calculation  tines  that  can  be  cxcccted  are 
appropriate.  The  tines  given  here  are  the  average  tines  used  in  com- 
puting the  results  given  above  and  ore  based  on  a computer  code  run  on 
a CDC  6600  computer.  The  average  computation  tine  was  1.5  sec  for  a 
5 shake  calculation  of  the  electric  field  development  computed  every 
.1  shake  using  50  steps  in  the  numerical  integration.  Increasing 
either  the  parameter  of  50  steps  or  5 shakes  by  any  factor  increased 
the  computation  tine  by  the  sane  factor.  The  presence  of  a precursor 
burst  had  no  noticeable  effect  on  computation  tine.  When  the  approxi- 
mations of  Eqs  (61)  and  (62)  were  valid,  the  tine  dropped  to  less  than 
.2  sec  for  a 5 shake  calculation. 

An  interesting  observation  related  to  the  numerical  integration  is 
that  as  Eqs  (61)  and  (62)  become  valid,  the  numerical  integration  could 
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V . Conclusions  and  RecovroncJo t i ons 

This  report  has  presented  a ealeu lot ional  model  for  the  hioh  al- 
titude EtfP.  This  model  falls  between  the  very  simple  and  niqhly  sophis- 
ticated models  currently  available.  The  fnodel  is  easy  to  understand 
but  allows  the  variation  of  the  major  theoretical  parameters.  The  model 
is  also  simple  enough  to  use  so  that  a preionization  level  car  be  in- 
corporated with  a very  minimum  of  effort.  The  results  given  by  this 
model  are  in  excellent  agreement  with  the  available  data  from  the  very 
sophisticated  models  usina  the  same  geometry.  The  computation  tine  is 
on  the  order  of  a few  seconds. 

Uses 

There  are  two  primary  uses  seen  for  this  model.  The  first  Is  for 
sensitivity  determinations  for  a given  parameter,  such  as  precursor 
yield  in  preionization  studies,  height  of  burst  for  a fixed  yield,  or 
changes  in  the  magnitude  and  inclination  of  the  neonagnetic  field. 
Calculations  using  the  HAEMP  model  can  help  to  evaluate  the  relative 
importance  of  the  many  variables  encountered  in  EMP  problems.  It  is 
useful  in  this  role  because  it  oives  quick  answers  which  arc  reasonably 
accurate. 

The  second  use  is  as  a classroom  model.  The  basic  physics  of  the 
high  altitude  problem  is  present  in  an  easy  to  use  form  which  demon- 
strates how  the  currents,  conductivity,  and  electric  field  interact 
with  each  other.  Because  the  calc  si  t ions  reauired  to  get  an  answer  are 
minimal  and  straightforward,  the  physical  processes  are  not  lost  in  the 
search  for  a result,  as  can  happen  all  too  easily  with  hiqhly  sophis- 
ticated models.  Since  the  results  are  in  relatively  good  agreement  with 
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the  sophisticated  models,  the  student  can  leave  the  classroom  with  a 


oood  concept  of  the  overall  problem. 


Limitations 


The  model  does  have  its  limitations,  however.  These  are  due  to  the 
various  assumptions  and  approximations  used  to  arrive  at  the  final 
calculat ional  model.  These  assumptions  and  approximations  should  be 
kept  in  mind  when  using  the  model.  If  this  is  done,  the  results  from 
any  problem  can  be  reasonably  interpreted. 

Two  of  the  error  sources  are  the  small  angle  approximations  and  the 
assumption  of  a flat  earth.  For  the  small  angle  case  the  error  Is 
nealigible  at  very  early  times  but  increases  as  the  time  increases. 

In  the  flat  earth  assumption,  any  error  is  relatively  constant  for  a 
particular  geometry,  but  changes  when  the  geometry  changes.  This  error 
increases  as  the  ground  inane  of  the  target  moves  on  the  earth's  sur- 
face away  from  the  ground  zero  ooint.  The  other  assumptions  also  add 
some  error  to  the  results,  but  if  the  effects  of  these  assumptions  are 
renembered,  the  KAEMP  model  can  be  helpful  for  the  solutions  of  a wide 
range  of  EHP  problems. 


P.ecomnenda  1 1 ons 

There  are  several  areas  of  investigation  which  might  lead  to  Improve- 
ments in  the  HAEMP  model.  One  of  the  more  important  possible  investi- 
gations is  to  determine  if  the  electron  transport  equations  used  In  the 
AFWL  CilEHP  model  to  describe  the  secondary  electrons  can  be  greatly 
simplified  and  used  In  the  HAEIIP  model.  Another  area  is  the  possibility 
of  allowing  recombination  and  cascading  of  secondary  electrons.  The 
use  of  self-consistent  electromagnetic  fields  Is  also  a possible  source 


1(1 


of  Improvement.  This  change  however  would  have  to  he  examined  very 
closely  to  determine  the  effects  it  would  have  on  the  validity  of  the 
small  anqle  approximations,  which  ore  vital  to  the  IIAEMP  model.  An 
additional  possibility  is  the  consideration  of  small  anole  scattering 
of  the  Compton  electrons.  This  could  possibly  be  included  by  a random 
elimination  of  some  Compton  electrons  for  Compton  current  production 
while  still  allowing  all  Compton  electrons  to  produce  secondary  elec- 
trons. 

Each  possible  change  to  the  HAEHP  model  should  be  examined  with 
consideration  to  the  tradeoffs  between  more  complete  solutions  and  the 
continued  ease  cf  use  and  understanding  of  the  model. 
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HAplP^odp  I Stun  nary 


The  user  of  the  HA£MP  mode!  must  first  determine  c he  following 


parameters  in  MKS  units: 


the  gamma  yield  Y of  the  weapon  in  MeV. 


the  height  of  burst  (hOB)  of  the  weapon. 


the  magnitude  30  and  angle  of  inclination  of  the 


geomagnetic  field  in  the  absorption  region. 


the  target  location  in  spherical  coordinates  with  the 


origin  at  the  burst  point  and  wi th  trc  polar  axis 


parallel  to  the  geomagnetic  field  lines. 


the  kinetic  energy  Ec  of  the  Comotor,  electrons  In  MeV. 


the  parameters  defining  the  gamma  output  pulse  shape. 


the  atmospheric  scale  height  S to  be  used. 


From  these  parameters,  the  values  for  the  annle  A and  the  two  dis- 


tances ®?1IM  and  RMAX  used  for  integration  over  the  absorption  reoion 


can  be  found  from  geometry.  The  angle  from  the  vertical  to  the  line  of 


sight  from  the  burst  point  to  the  target  is  the  angle  A,  as  shown  in 


-ig.  1 on  page  5.  The  value  of  RHIM  is  the  radial  distance  from  the 


ourst  point  to  the  top  of  the  absorption  layer  (50  kn  above  the  earth’s 


surface)  along  the  line  of  sight  to  the  target.  The  value  of  RMAX  is 


the  corresponding  distance  to  the  bottom  of  the  absorption  layer  (20 


Km  above  the  earth’s  surface).  See  Fig.  1 on  pane  5. 


The  electric  field  development  from  tines  0 to  t in  steps  of  At  in 


the  retarded  tine  frame  is  found  by  the  numerical  integration  of 


1 7 If  (rEe>  + uo  J8  + "o  o(t)  e8  " 0 
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"fir  (rE0>  + "c  \ + »c  «<1>  E„  • 0 (* 3) 

over  r from  RMIM  to  RMAX  for  each  tire  increment.  The  initial  condi- 
tions on  these  equations  are  that  E9  a rd  E?  are  zero  at  t = 0.  The 
number  of  steps  taken  over  the  inteeration  interval  will  depend  on  the 
yield  of  the  weaoon  and  the  accuracy  desired,  ho  less  than  50  steps 
should  be  used  or  more  than  500  steps  necessary  for  gamma  yields  up 
to  100  Kt. 

The  solutions  to  Eqs  (k2)  and  (A3 ) will  give  the  electric  field 
components  at  r <*  RMAX.  The  total  field  strenoth  Is  oiven  by 

E = [(E0)2  + (Eg)2]^  (63) 

and  the  field  strength  at  the  target  is  given  by 

, RMAX  E 

E (target)  « (53) 

'target 

The  pe.iK  electric  field  value  over  the  time  interval  0 to  r Is  the 
maximum  value  found  using  Eq  (53). 

In  the  special  case  where  the  conditions 

If  •!?  (rEeyi  « “0  iJe,jii  (59) 

cf  if?  (rE8,e>l  c<  b0  o(t)  E9>0  (60) 

hold,  the  electric  field  corcoonents  are  alven  by 

Efl  - - Jfl.Mt)  (61) 


Eg  ‘ • Jjj/d(*) 
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fhc  values  from  tqs  (61)  and  (6?)  nay  then  be  used  with  Eqs  (6A) 
ar.d  (53)  to  find  the  electric  field  at  the  toroet. 

In  order  to  integrate  tqs  (A?)  ar.d  (C3)  or  to  compute  tqs  (61) 
and  (62),  the  terms  J§,  J§,  and  o(t)  are  needed  as  functions  of  time  i 
and  distance  r.  These  are  given  by 


j!i  = eq(r)  Vc  sin  0 cos  8 — j t'j  f(T)  dt' 


Jg  = - og(r)  VQ  sin  8 a | t'  f(T)  dt' 
**  o 


o(t)  = n^nS<T> 


Y 1 s / HOB  x - rrcosAv  ,nl 

c exo  . - T — exP  * ' ~ T*  , exr,  , c , * 1 L 

EC  W2  A (r)  I.  x0  COS  A ^ S / >_  V S J }i 


f(T)  = [exp  <-aT)  - exo  (-bT)]  U(T) 

b - a 


QV_  t R/V 

n«(t)  ■ — — g(>‘)  j [j  f(T')  dt"]  dT- 


t = t - (i-e)  t' 


T'  = t'  - (1-S)  t' 


withe  - ■£.  ar.d  q = Ec/33  eV 


The  Compton  electron  cyclotron  frequency  u is  given  by 


1|7 
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r 


F: 


! 

i 


i 


1 


L 


The  Conoton  electron  velocity  V0  is  oiven  by 


V 


o 


c 


( — !V 

^Ec  + m0  > J 


The  Comotoo  electron  range  at  sea  level  is  given  by 


(53) 


R = (612  £c<’-265  - ’W'  Ec)  )/D  (5, , 

where  9 = 1.293  mq/cm3.  This  range  is  then  converted  to  meters  and 
scaled  to  the  appropriate  value  for  the  point  r by  multiplying  the  com- 
puted value  by 

exp  {(HOB  - r cos  A)/S]  (66) 

Then  the  Compton  lifetime  P/V0  is  oiven  by  the  quotient  of  the  correc- 
ted P and  V0,  but  is  never  allowed  to  exceed  one  usee. 

The  value  of  A(r)  is  given  by  the  product  of  (300  m)  and  Eq 

(63). 

The  collision  frequency  vc  at  sea  level  is  given  by  the  maximum 
of 


vc  = -2.5  do)20  t + 6. 5 do)12 

(55a) 

vc  = 0.63  (i0)8  E + 1.6  (10)12 

E < 5 do)1- 

(55b) 

vc  = 6 (10)7  E + 0.8  (10)12 

E > 5 (10)“ 

(55c) 

or  2.3  (10)12  where  t is  the  retarded  time  in  seconds  and  E is  the  elec- 
tric field  strength  in  volts  per  meter.  The  collision  frequency  is  not 
allowed  to  be  greater  than  6.6  (1 0) 12 . The  value  of  the  collision  fre- 
quency is  recomputed  after  each  tine  step  in  the  development  of  the 

<l3 


,J 


I 


gep/pi</75-13 


electric  field.  The  sea  level  value  is  also  scaled  for  altitude  by 
dividin'}  by  Eq  (of* ) . 

The  intcqrals  in  Eqs  (24),  (25).  and  (19)  can  be  evaluated  in 
closed  form.  The  steo  function  II (T)  in  Eq  (43)  results  In  two  expres- 
sions for  each  of  the  inteorals.  The  first  set  of  exnrcssions  Is 
used  when  r — (?  - 2)  (R/V0)  <.0  and  the  second  when  x * (I  - S)  (R/VQ) 

> 0. 
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Jc(t)  =>  eo(r)  V sin  6 cos  0 


(I  - K)3 


{exp  (-at)  {exp  (aL)  [(aL)2  - 2aL  + 2)  - 2) 


exp  (-bt)  (exp  (bL)  [ (bL) 2 - 2bL  + 2]  ~r) 

b* 


J§(t)  = - cp(r)  V0  sin  6 w 


b - a (,  . B)2 


{exp  (-at)  {exp  (aL)  [aL  - l)  + i] 


- exo  (-bt)  (exp  (bL)  f bL  - 1)  + 1]  r) 


o(t)  = -~g(r)  -rrn,rs 


(|-  [at  - 1 + exp  (-aL)  - [exp  (aL)  - l][exp  (-at)  - exp  (-aL)]] 


- bL  - 1 + exp  (-bL)  - [exo  (bL)  - l][exp  (-bT) 
b 


- exo 


where  L * (1  - 3)(R/V0). 


50 


GLP/PH//5-13 


^ 1 


Vita 

Louis  W.  Seiler,  Jr.  was  Lorn  on  1C  January  I9L7  in  ilew  Orleans, 

LA-  He  araduated  from  high  school  in  Ilew  Orleans  In  196<i.  He  attended 
Louisiana  State  University  in  Kew  Orleans  and  received  the  deqree 
Bachelor  of  Science.  Physics  in  19f>S-  He  entered  the  U.S.  Air  Force 
Officer  Training  School  in  October  19C9  and  received  a commission  in 
December  1969-  He  served  as  Group  Interactive  Computer  Systems  Analyst 
in  the  fiOOO  Aerosoace  Applications  Group  until  June  1973  when  he  entered 
the  Air  Force  Institute  of  Technology. 


Permanent  Address:  35 1 S Crovdon  Street 

Slidell,  LA  70h5? 


This  thesis  was  typed  by  Katherine  Randall 


51 


